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We study the linear Hall response of 2D ballistic system on inhomogeneous magnetic field. We 
establish that in classical limit the Hall conductivity response on local magnetic field is quantized 
in units of an = , lie ■ The quantized value depends on the point where the field is applied, this 
dependence being irregular in chaotic billiards. The phenomenon allows for direct tracing of special 
electron trajectories that belong to fractal repellor. We discuss how quantum effects smooth the 
quantization. 
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' Since the first ballistic constriction in 2D electron gas has been formed |l[ a great deal of experimental and theoretical 
attention has been devoted to magnetotransport in such microstructures (see ^ for review). Several effects such as 
mesoscopic fluctuations can only be understood on quantummechanical basis. However the main features of 
magnetotransport are to be comphehended in classical framework j^] and can be naturally explaned in language of 
classical electron trajectories. 

c/3 [ An interesting new development is to use the ballistic constriction as a sensitive detector of inhomogenious magnetic 
^ • field. 01 In this experimetal work, the constriction has been used to reveal the details of magnetization of a microscopic 
G ' superconducting particle. However, the intrinsic physics of this detector also appears to be interesting and complicated. 
There are several theoretical and numerical works tackling the problem of inhomogenous magnetic field. |^] In the 
present paper we adopt the classical approach and concentrate on the Hall response linear in inhomogeneous magnetic 
field applied. 

'■^ , We establish two main results: i. the linear Hall response is quantized in units of an = 2^^hc ^ quantized value 
G ' depending on the position in the constriction ii. The lines that separate regions of different quantized values outline 
O , the electron trajectories that belong to fractal repellor. Q The fractal repellor is a manyfold of electron trajectories 
orbiting inside the constriction for infinite time. This fact is of importance since it opens the way for direct observation 
of the chaotic nature of scattering in a ballistic constriction in a transport measurement. This is hard to believe since 
any transport quantity results from averaging over extended regions of phase space. Thus one could expect that the 
chaotic dependence of the scattering result on initial conditions is averaged out. 

Let us start with the theory formulation for arbitrary multi-lead ballistic constrictions. Q (Fig. 1). The currents 
Ii in the leads are related to the voltages by means of conductivity matrix: Ii = dj Vj . It obeys Onzager relations 
C Gij (H) = Gji{—H) and provides conservation of current: ^ ^ Gij = 0. In classical limit, the conductivity is determined 

from the master equation 



df p df dU(r) df e , 9/ „ 

' ^ ^ ' -(p X H) • = (1) 
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S for distribution function /(r,p). We concentrate on deviations of the conductivity matrix which are linear in H. 
■ Quite generally, those can be presented if the form: 

I 2 « 

3 : SG,, = J d\N,,{v)H{v). (2) 

• , Now we are to demonstrate that Nij can only take integer values. We expand distribution function in the form 
/ = ./o + /ii /i being proportional to the magnetic field, that gives equation for /i: 

at m or or Op mc Op 

We can find the current generated by /i without explicitly solving Eq. ||. To do this, we note that the right hand 
side of (1) can be presented as a sum over particle sources. We introduce the probability Pi{x, y, p) to get to the lead 
i while starting in point x,y with momentum p. Then the current deviation can be presented in the form 
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Now we shall specify /q. Since the system is ballistic, /o is contributed by electrons coming from the leads. Since the 
probability to come from i — th lead is Pi{x, y, — p), fo{x, y, p) = X]j Pii^i ~P)/F(e(p — eVi). The deviation of the 
current is thus given by 

"•-^S / J^«(^.P)(P-)^^S|_P);H.(p-.va)//M (5) 

We note that the electrons participating in transport are concentrated in a narrow energy strip near Fermi level. 
That allows us to disregard the dependence of P upon energy or modulus of p. P becomes now a function of angle 
(f) that determines the bearing of p. Under this assumption, Eq. ^ can be integrated over energy. We use the density 
of states V = m/nh^. This finally yeilds Eq. ||with 



where the angle parametrizes direction of p. Since the motion of electrons is classical, P{4>) can only take two values: 
and 1. So that Nij is integer. We stress that this proof uses only the fact that the scattering in the constriction is 
a classical one. It does not depend on details of potential relief. Also, the minor modification of the derivation shows 
that the result holds for an arbitrary electron dispersion law provided the Fermi surface remains single-connected. 

Does this mean that the linear Hall response is physically quantized? No, not necessarily. The physical, observable 
quantization implies that N{x,y) not only takes integer values at any point but also retains the same value in some 
finite neighborhood of this point. We will see below that sometimes it is not the case. 

Let us have a closer look at functions Pi{(j)) , Pj {(j) — tt) (Fig. 1). We consider, without loss of generality, only 
constrictions with smooth edges. Then Pi((/)) develops a typical self-similar chaotic stucture. Q Although it can only 
take two values, and 1, it changes infinite number of times between these values. Moreover, in any finite neighborhood 
of each change one finds infinitely many other changes. In order to quantify the situation, one subdivides the whole 
interval of angle on smaller intervals of the length 2tt/K, K being big integer number. An interval can be "black" 
if Pi = 1 for all its points, "white" if = O for all its points, or "gray" if there is at least one change within the 
interval. With increasing resolution K, the total number of "gray" intervals M scales as 

M ~ K"^" , (7) 

dn being Hausdorf dimension of the fractal manyfold of changes. The dimension dn does not depend on a point 
where Pi is measured. It is determined by the overall geometry of the constrictions. 

Let us now analyze the expression at finite resolution K. First it seems that the answer is not determined at 
any finite K since the answer may depend on the behavior of P within "gray" intervals. However, we notice that 
the expression is completely determined provided all "gray" intervals of Pi are covered by "black" or "white" 
intervals of P, and, similar, all "gray" intervals of Pj are covered by "black" or "white" intervals of Pi. This stems 
from the structure of integrand: if Pi is a constant within an interval (a, 6), the contribution of this interval equals 
Pi{Pj{a) — Pj{b)) and docs not depend on the details of Pj. Let us assume that for well-developed chaos positions of 
"gray" intervals for Pi and Pj are statistically independent. We estimate the probability for all "gray" intervals be 
covered. Since the probability to find a "gray" interval against a given interval is M/K, and the number of "gray" 
intervals is M, we obtain 

Pco.er = {l-^r ^OM-^) (8) 

in the limit of M oo, M/K 0. Since M ~ K"^" , we find two distinct regimes. If dn > 1/2 Pcover approaches 
with increasing resolution. This means that N(x, y) can not be an integer in any finite region of x, y but rather shall 
exhibit continuous behavior. In opposite case dn < 1/2 Pcover reaches 1 so that N{x,y) takes only integer values 
that persists in finite regions. 

Let us now reformulate the above statements stressing the relation between the changes of Pi,j{x,y, (f>) and the 
trajectories of electrons that exhibit chaotic scattering in the constriction. Q The trajectory that starts at the point 
(a;, y) with the bearing (j) at which a change occur may end neither in the lead i (Pi = 1) nor in any other lead {Pi = 0). 
Rather, it stays in the constriction for infinite time. Let us now consider the trajectory of opposite bearing, (f) + n. If 
for all infinite trajectories their counterparts of opposite bearing are finite, we recover the above mentioned situation 
where all "grey" intervals are covered by "white" or "black" intervals. Thus N{x,y) is defined and integer in this 
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point. Albeit this is not the case if in a given {x^y) at some (j) there is a trajectory that is infinite for both bearings, 
(j) and TT + (/>, so that N{x^ y) remains undefined in this point. Owing to time reversal symmetry, this trajectory must 
correspond to a periodic orbit in the constriction. These periodic orbits form what is called a fractal repellor of the 
constriction. Therefore the regions of constant N{x,y) can be separated from each other only by the trajectories of 
the fractal repellor. This is how these trajectories can be seen. 

An imaginable visualization experiment would resemble those with STM or AFM. Let us bias the constriction and 
measure its Hall response. As a source of magnetic field we can take a long and sharp magnetic needle with one of 
its ends hanging close to the sample surface. We can scan the whole constriction area with this needle measuring 
the Hall response as a function of (x, y) needle position. If we contrast the changes of the response we will see 
bunches of the fractal repellor trajectories. (Figs. 2-4) 

To illustrate the visualization of the fractal repellor, we simulate classical electron trajectories in the cross-like 
billiard with smooth edges (insert in Fig. 2), these edges being quarter circles. This system resembles the actual 
experimental setup. Q If the Hall conductivity is defined in a conventional way as the ratio of Hall current to the 
voltage drop across the whole constriction, it is quantized in units of an /2. The chaotic properties can be characterized 
by a single dimensionless parameter /3 = R/d, R being circle radius, d being the lead width. This parameter can 
be in principle changed in situ by means of gate electrodes. We stress that such a system corresponds to a classical 
model problem of chaotic scattering: scattering by unpenetrable discs. Although a three-disk billiard received most 
of the attention Q , the four-disk billiards have been thoroughly studied as well. [D The Hausdorf dimension denined 
as above was shown to change monotonically between 1/3 {/3 = 0) and 1 (/3 ^ 1) with increasing f]. The transition 
between discrete and continuous behavior of N{x, y) takes place at /3 w 0.5. We have chosen /3 — 1/3 that corresponds 
to dH « 0.4. 

We have calculated approx. 100000 periodic orbits, some of them being equivalent due to symmetry of the problem. 
We plot them in Fig. 2, 3, 4 zooming on the constriction center with zoom factor of 10. We observe self-similar 
pattern which is characteristic for chaos. Despite the huge number of the trajectories, and dn close to 1/2, they do 
not uniformly cover the whole area of the constriction. Rather, they are grouped in bunches with well resolved space 
in between. We plot the local Hall conductivity along y symmetry axis in right plane of each figure. That allows one 
to see the correspondence between trajectories of the repellor and quantized values of N{x, y). 

Quantum mechanics is known to smooth peculiarities of the classical chaos. |8| Although the representative quantum 
mechanical solution of the problem is hardly within the reach due to big number of basis states required, we can draw 
simple estimations of quantum effects in the billiards. We assume that the angle resolution K/2n is limited by a 
typical diffraction angle — ^JX/L, Here A, L are the electron wavelength and the typical size of the constriction, 
respectively. This implies that the regions of quantized values of N{x, y) that would require better angular resolution, 
remain unresolved and N{x^ y) is continuous there. Eq. ^ gives an estimation of relative area of these unresolved 
regions, 

A5/5~exp(-(L/A)i/2-''«). (9) 

In conclusion, we have predicted that the local Hall conductivity of classical ballistic constriction may be quantized 
in values of an = j^^rj—- This may open up new ways to measure fundamental constants. There is a transition 
between quantized and continuous behavior of the local conductivity that takes place at dn = 1/2. The regions of 
the same quantized N{x,y) are separated by trajectories of the fractal repellor, that provides the means for their 
visualization. 
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[1] B. J. van Wees et al. Phys. Rev. Lett. 60, 848 (1988). 

[2] C. J. W. Beenakker and H. van Houten,in:i?teciromc properties of multilayers and low-dimensional semiconductor structures, 

Ed. by J. M. Chamberlain et al.. Plenum Press, New York (1990). 
[3] B. L. Altshuler, Pis'ma Zh. Eksp. Teor. Fiz. 41, 530 (1985) [JETP Lett. 41, 648 (1985)]; P. A. Lee and A. D. Stone, Phys. 

Rev. Lett. 35, 1622 (1985); B. L. Altshuler and D. E. Khmelnitskii, JETP Lett. 42, 559 (1985). 



3 



[4] A. K. Gcim ct al., Suporlatticos and Microstmctures, 23, 151 (1998); A. K. Geim et al., Appl. Phys. Lett. 71, 2379 (1997). 

[5] X. Q. Li, F. M. Pootors, and A. K. Gcim, Journ. of Phys. C 9, 8065 (1997). 

[6] B. Eckhardt, J. Phys. A: Math. Gen. 20, 5971 (1987). 

[7] M. Buttiker, Phys. Rev. Lett. 57, 1761 (1986). 

[8] P. Gaspard and S. A. Rice, J. Chem. Phys. 90, 2225 (1989). 

[9] P. Gaspard and D. Alonso, Phys. Rev. A 45, 8383 (1992); P. Gaspard, D. Alonso, T. Okuda, and K. Nakamura, Phys. Rev. 
E 50, 2591 (1994). 

FIG. 1. Multi-lead ballistic constriction. The Hall response on magnetic field in the point {x, y) is determined by probabilities 
Pij{x,y,<j))- Their angular dependence is sketched on the right. Grey intervals comprise fine chaotic structure of P that can 
not be resolved at a given angular resolution. 

FIG. 2. The fractal repellor trajectories in the constriction. Only a quarter of the constriction is shown. The right panel 
presents quantized linear Hall response along the y symmetry axis. 



FIG. 3. The same as in Fig 2. but at 10 times smaller scale. 



FIG. 4. The same as in Fig 2 at 100 times smaller scale. 
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